Reducing Disorder in Artificial Kagome Ice 
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Artificial spin ice has become a valuable tool for understanding magnetic interactions on a mi- 
croscopic level. The strength in the approach lies in the ability of a synthetic array of nanoscale 
magnets to mimic crystalline materials, composed of atomic magnetic moments. Unfortunately, 
these nanoscale magnets, patterned from metal alloys, can show substantial variation in relevant 
quantities such as coercive field, with deviations up to 16%. By carefully studying the reversal 
process of artificial kagome ice, we can directly measure the distribution of coercivities, and by 
switching from disconnected islands to a connected structure, we find that the coercivity distribu- 
tion can achieve a deviation of only 3.3%. These narrow deviations should allow the observation of 
behavior that mimics canonical spin-ice materials more closely. 



Water ice and spin ice are classic examples of geometri- 
cally frustrated systems P, Q , both with residual low-T 
entropy In water ice, thermodynamic phases with 
orderedprotons were discovered after decades of experi- 
ments y. In contrast, no dipole-ordered phase has been 
observed in spin ice even at the lowest accessible temper- 
atures, contrary to a theoretical prediction (Hj. Divergent 
relaxation times and quenched disorder in samples have 
been cited as possible explanations. Artificial spin ice 
has been proposed to help address these questions (6| , as 
it allows the direct control of the geometry of the lattice, 
with the combined ability to directly image the result- 
ing microstate. Here, samples are composed of lattices 
of nanoscale ferromagnetic islands, where the magnetiza- 
tion of each element points along its longitudinal axis. At 
the vertices of the lattice, the ferromagnetic elements in- 
teract, and because of the geometry of the system, their 
magnetic configurations are frustrated This al- 

lows the study of frustration in systems where crystalline 
imperfections can be completely removed by design, or 
introduced in a controlled way. Unfortunately, current 
lithographic techniques are limited by unintended rough- 
ness at edges and interfaces, creating inadvertent disor- 
der. This diminishes the ability to compare observations 
from artificial spin ice materials with studies of spin ice 
oxides, where magnetic atoms are presumed to be iden- 
tical. 

Edge roughness of nanomagnetic elements is known to 
substantially influence the coercive field, by creating nu- 
cleation sites that can initiate the magnetic reversal . 
In some artificial spin ice geometries, this edge roughness 
can create a large variability in the behavior of the arti- 
ficial "atoms" (magnetic nano- islands ) . In recent studies 
of artificial kagome ice, the variations in coercivity were 
found to be substantial — up to 16% of the average co- 
ercive value d, 0, [13] . This variability can easily be re- 
duced by choosing materials with low crystal anisotropy 
0, but we here show substantial further reduction with 
a geometry with connected magnetic islands. In a con- 
nected geometry, nontrivial spin textures (domain walls) 



already exist at the vertices even in equilibrium, and thus 
the reversal process is not dictated by the nucleation of 
new domain walls, a step that is more susceptible to the 
effects of edge roughness. 

To address the problem of disorder and variability in 
artificial ice, techniques must first be established to mea- 
sure the static disorder in the material, a problem that is 
complicated by the sought-after occurrence of statistical 
disorder, in analogy with thermal disorder in pyrochlore 
spin ice. Prior work has addressed this problem by focus- 
ing on the reversal of artificial kagome ice in an external 
field applied at 180° to the direction of initial magnetiza- 
tion [Sl, |9|, a potentially deterministic process, in which 
statistical variations are hopefully minimized. However, 
the reversal process under these conditions shows large 
discrete avalanche events, and for the materials with 
small static disorder described below, the variations in 
these reversals are almost entirely dominated by non- 
reproducible statistical disorder. Prior studies provided 
disorder estimates as large as 16% through Monte Carlo 
modeling of the 180° reversal process, but this approach 
is not possible when the amount of quenched disorder is 
much smaller. 

In this Letter, we present a method to alleviate this is- 
sue by performing reversals with the field applied at 120° 
and 100° to the direction of initial magnetization instead 
of 180? Under these conditions, magnetization reversal 
in individual islands proceeds largely independently of 
their neighbors, preventing the formation of avalanches. 
This enables us to extract the spread of coercivities di- 
rectly from the magnetization curves M{H). These data 
are then used to calculate a disorder parameter, a/Hc^ 
where He is the average coercive field of a magnetic ele- 
ment and fj is the standard deviation of the distribution. 
Using this parameter, we are able to directly compare our 
data with those from other groups. By considering care- 
fully the magnetic reversal process, we demonstrate that 
this disorder parameter can be greatly reduced in a sys- 
tem with connected magnetic islands. Additionally, all 
data used to determine the disorder present in a crystal 
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can be obtained directly from experiments rather than 
relying on parameters extracted from Monte Carlo sim- 
ulations. 

In our studies, we choose the kagome lattice for our 



artificial spin ice [sMllj. as shown in Fig. [UJa) and[ljb). 



Samples were patterned via e-beam lithography onto an 
electron transparent SiN membrane with PMMA resist 
llj . The Ni8oFe2o films were deposited by e-beam evap- 
oration followed by metal lift-off. The resulting elements 
are 500 nm long, 110 nm wide, and 23 nm thick with 
an additional 4 nm capping layer of Al to prevent oxida- 
tion. Reversals were performed in situ in a transmission 
electron microscope in Lorentz mode to capture the mag- 
netic configuration of the crystal throughout a reversal. 
A magnified image of a portion of the sample can be 
seen in Fig. [Ija). Magnetic reversal in this system can 
be described in terms of the emission, propagation, and 
absorption of domain walls carrying non-zero magnetic 
charge 1^. Similar emergent magnetic monopole exci- 



tations have been shown to exist in both conventional 
iB - 20l | and artificial spin ice [1, 0] . 

The nature of magnetic charges is understood as fol- 
lows. Because magnetic induction B = ^o(H -I- M) is 
divergence-free, a region where lines of magnetization M 
terminate or originate becomes a source of sink of mag- 
netic field H. Thus magnetic charges in artificial spin 
ice contain integer multiples of the magnetization flux of 
a single magnetic element. In these units, a monopole 
propagating along an element during reversal can carry 
a charge ±2, and a vertex can carry a charge of ±1 or 
±3. The interactions between the charges via H can be 
described simply in terms of a Coulomb potential fl3 |. 
Triple charges represent ice- rule violations for the kagome 
lattice, and although they have been observed experimen- 
tally during 180° reversals in other systems |8), iB], they 
never appear in the low-disorder system we study here 

Before considering the angle dependence, it is instruc- 
tive to consider the microscopic details of the reversal 
process, depicted in Figure [TJc). In panel (z), the hori- 
zontal element will be the first to reverse because the ex- 
ternal magnetic field has the largest component along its 
axis. The reversal begins with the emission of a domain 
wall of +2 charge into the horizontal element, leaving be- 
hind a vertex with —1 charge. The domain wall moves 
along the element until it reaches the opposite vertex, 
thereby reversing the element's magnetization. In panel 
{Hi), the vertex with charge —1 pulls the domain wall 
with charge +2 in, thereby changing its own charge to 
-|-1. The reversal is different for the disconnected lattice 
where a domain wall needs to be nucleated at the end of 
an element, rather than emitted from a vertex. In this 
case, the field required to inject a charged domain wall 
can vary significantly depending on crystal symmetry and 
edge roughness of the elements. 

One might expect the magnitude of the coercive field 




FIG. 1. (a) In-focus TEM image of the artificial kagome ice. 
Lines representing the kagome planes are overlaid to show 
how the origin of the honeycomb structure of our crystals, (b) 
Lorentz contrast image of the kagome lattice. The contrast 
difference can clearly be seen with a bright and dark edge 
across each element in the array. The magnetization of six 
elements is shown along with the resulting charge at the vertex 
they are connected to. Inset: the intensity profile along the 
indicated black line. The magnetization of two elements is 
captured in the profile and depicted by arrows on the Lorentz 
image and on the profile itself, (c) A cartoon of the reversal 
process of an individual element is shown in panels (i)-(ro), 
as described in the text. 
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FIG. 2. Illustration of the crystal magnetization at the be- 
ginning oi sl 8 — 120° reversal. The table shows values of 8i 
for the elements of the three sublattices. The angles are with 
respect to Happ- The red lines at the vertices indicate the 
location of a head-to-head or tail-to-tail domain wall which is 
the source of magnetic charge at each vertex. 



to be lowest when the field is parallel to the element 
receiving the emitted domain wall: only the longitudinal 
component of the magnetic field pushes a domain wall 
along the element. However, in a connected lattice, the 
transverse field component also plays a role, thanks to an 
asymmetric distribution of magnetic charge within each 
vertex. 
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Consider a simplified model of a vertex consisting of 
three thin domain walls meeting at the center of the ver- 
tex, Fig. [5] In a vertex with unit charge, two of the walls 
are neutral (head-to-tail) and one is magnetically charged 
(head-to-head or tail-to-tail). An applied field exerts a 
maximum force on the charged domain wall when it is 
applied normal to the wall, or at 30° to the axis of the 
element receiving the wall. Micromagnetic simulations 
confirm these qualitative considerations The coer- 

cive field depends on its orientation as follows: 



H,i0,) = WJ\cos{ 



(1) 



Here 6i is the angle of the applied field with the long axis 
of the z*'' element, is the minimum coercivity of that 
element, and a = 19° is an offset angle [l^ . 

Using this model, we measure the disorder in our arti- 
ficial ice, as caused by variability in the coercive field of 
the elements, by performing magnetic reversals of large 
arrays with more than 10* elements. The system is ini- 
tially magnetized in a high field H ^ He- The field is 
switched off and the sample is rotated through the de- 
sired angle 6 in its plane. In this Letter, we focus on 
e = 100° and 120? As can be seen from Eq. (P), these 
orientations guarantee large differences in the coercive 
fields for different sublattices and lead to avalanche-free 
reversal. This allows us to extract reversal statistics for 
both sublattices during the reversal process. 

During reversal, the field is increased in steps of 5.0 
Oe, and since it can also effect the Lorentz imaging, the 
field must be reduced to after each step, to capture 
an image at remanence. Images are then run through a 
MATLAB script that automatically determines the mag- 
netic configuration of each element; one image will give 
us one point on the M{H) plot, as shown in the top panel 
of Fig. [31 These results can then be compared with nu- 
merical simulations of an array of interacting elements as 
described previously Simulated curves M{H) for a 
Gaussian distribution of coercive fields with a /He = 5% 
are shown in the bottom panel of Fig. [3l D etails of the 
simulations will be described elsewhere |2l| . 

To determine the disorder in our crystals directly, we fit 
the M{H) curves to the expected cumulative distribution 
function, a superposition of error functions in the case 
when coercive fields have Gaussian statistics: 



M4H) 



(2) 



where Mx is the magnetization along the initial polar- 
ization direction, i indexes the sublattices in order of re- 
versal, rrii is the amount of magnetic moment carried 
by the sublattice, is the average coercive field, and 
ai is the standard deviation. From fitting the reversal at 
6 = 100° we obtain a disorder parameter, (Ji/Hl = 3.3%, 
from the values Hi = 347 Oe and ai = 11.3 Oe. 
The second plateau at 6* = 100° gives a slightly higher 
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FIG. 3. Experimental (top) and theoretical (bottom) results 
of a uniaxial reversal along 100° 120° and 180° Solid lines 
in the experimental results are fits of a superposition of er- 
ror functions used to extract the reversal statistics for each 
sublattice. 



disorder parameter of 4.7%, from H"^ — 915 Oe and 
(72 = 43.3 Oe. We believe this is due to a rotational disor- 
der from element-to-element in the lattice, as expressed 
by the a offset parameter from Equation Disorder 
in a should not be expected to contribute to Oi when 
|cos(6'i -f a)| w 1, which occurs for cti near 9 — 100° in 
our system. When 6i « ±90° , the coercivity diverges and 
even small rotational disorder would be expected to cause 
substantial variations. From the reversal ai — 120° de- 
grees, we obtain (Ti/HI — 3.7%, and (T2/H'^ — 4.6%, 
confirming the trend of minimum disorder near the mini- 
mum H^, and increasing monotonically with 9 away from 
that point. We note that the reversals ai = 180° have 
a different character. They involve magnetic avalanches 
in which long chains of links reverse together. We plan 
to address this question in a separate publication. 

From this, we find that the static disorder is very low 
in our samples, 3.3%, which is a factor of 4-5 smaller 
than values for other systems [1, This low amount 
of disorder is a result of the geometry of our connected 
elements. Here, charged domain walls always exist at 
the vertices, even at remanence. For a given magnetic 
configuration, there is very little change in the location 
and angle of the charged wall throughout the sub-lattice. 
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0.167 


Co, connected [8] 


Experiment 


0.130 


NiFe, disconnected [9| 




0.033 


NiFe, connected, this work 


OOMMF 


0.040 


NiFe, connected 


Simulations 


0.089 


NiFe, disconnected 



TABLE I. Experimental reversal results and their associated 
disorder. Theoretical results are shown for comparison. The- 
oretical results are averaged over 10 different simulations for 
each lattice. 



herein rigorously obeys its corresponding set of ice rules, 
a desirable feature given the interaction energy of ap- 
proximately 10"* K [6|. This presents the system as a 
more ideal model of frustration, providing the hope of 
addressing outstanding problems, such as the occurrence 
of a possible ground state in kagome ice 2^ 24 1. 
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The same is not true for a disconnected lattice. The 
process of nucleating and injecting a domain wall is more 
sensitive to quenched disorder in disconnected lattices, as 
the ±2 domain walls must be nucleated from the end of 
one of the elements instead of simply injected, and any 
significant edge roughness may influence the nucleation 
process. This leads to an increase in the width of the 
static disorder for a disconnected lattice. 

We summarize the effect of the connected lattice on 
the disorder in Table U where we compare a / He with 
two other experimental results, using different materials 
Q and disconnected geometry Q . The effect of the con- 
nected geometry is to decrease the disorder from 13% to 
3.3%. We acknowledge that the connected lattice pre- 
sented herein and the disconnected lattice used for com- 
parison have been fabricated by different research groups, 
but we have also performed simulations that support 
the relative difference. These micro-magnetic calcula- 
tions were performed using the OOMMF software pack- 
age from NIST 22|. Simulations were carried out for 
180° reversals on Y-junctions with random edge rough- 
ness, comparable to the fabricated structures, and the 
results from an ensemble of calculations are also shown 
in Table m exhibiting a similar relative relationship. 

In addition to the lattice geometry, the material used 
in the fabrication of the crystal is important. Our crys- 
tals and most other artificial spin ices have been made 
using Ni8oFe2o, which has a low crystal anisotropy, giv- 
ing polycrystal films with uniform magnetic properties. 
This is not the same for cobalt, as used in Q, which has 
a strong hexagonal anisotropy. Magnetic elements fab- 
ricated from polycrystal Co films would be expected to 
have a broad distribution of coercivities. 

By using a connected lattice of Ni8oFe2o elements, we 
have alleviated both of these problems and created crys- 
tals with low disorder. Our crystal geometry is ideal in 
that it allows for a more accurate representation of spin 
ice in frustrated magnetic materials. For example, unlike 
other systems [1, U, the low-disorder system described 
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